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ABSTRACT 

We report the discovery of a low-mass companion to HR 3549 , an AOV star surrounded by a debris 
disk with a warm excess detected by WISE at 22 p,m (lOcr significance). We imaged HR3549B in 
the L-band with NAOS-CONICA, the adaptive optics infrared camera of the Very Large Telescope, 
in January 2013 and confirmed its common proper motion in January 2015. The companion is at a 
projected separation of ~ 80 AU and position angle of ~ 157°, so it is orbiting well beyond the warm 
disk inner edge of r > 10 AU. Our age estimate for this system corresponds to a companion mass in 
the range 15-80 Mj, spanning the brown dwarf regime, and so HR 3549 B is another recent addition to 
the growing list of brown dwarf desert objects with extreme mass ratios. The simultaneous presence of 
a warm disk and a brown dwarf around HR 3549 provides interesting empirical constraints on models 
of the formation of substellar companions. 

Subject headings: stars: brown dwarfs, stars: low-mass, stars: imaging, instrumentation: adaptive 
optics, instrumentation: high angular resolution 


1. INTRODUCTION 

Although high contrast imaging of self-luminous exo¬ 
planets or brown dwarfs around young stars is already 
difficult because of the angular resolution and contrast 
required, a further complication is that the interpretation 
of the nature of any companions found in these images 
is made difficult because it relies upon theoretical mod¬ 
els of formation and evolution that depends on a good 
knowledge of the host star characteristics. While instru¬ 
ments, observing strategies, and post-processing tech¬ 
niques have undeniably improved over the past decade, 
accurate stellar age determination is becoming the most 
pressing challenge of characterization. Because they are 
a robust yet imprecise sign of youth, it is no coincidence 
that the majority of exoplanets imaged by adaptive op- 
tics are found in s ystems with debris disks: HR879 9 
(Marois et al.||2008|), 0 Pictori s (iLagrange et al.||2010|), 
HD95 U86 jRa ^au et ai.||2U13|), Fbmalh aut jKalas et al. 
20081), GjUb (Fl'ha lmann et al. 2009), and HDlObOOb 
( [Hailey et al.|[^14[ ). We note that the GJ504 system is, 
out of all bona hde directly imaged exoplanets, currently 
the only one without confirmed infrared excess, but is 
also the subject of an on-going controversy about its 


dmawet@astro.caltech.edu 

^ Jet Propulsion Laboratory, California Institute of Technol¬ 
ogy, 4800 Oak Grove Drive, Pasadena, CA 91109 

^ Based on observations made with ESO Telescopes at the La 
Silla Paranal Observatory under programs: 090.C-0486A and 
094.C-0406A. 


age. Fuhrmann & Chini (2015) indeed suggests that the 
GJ504 system is actually of Solar age, and that GJ504b 
is thus likely not a ~ 4 Jupiter mass planet, but a 4.5 
Gyr old Brown Dwarf. 

Debris disks are the signposts of planetary systems. 
Gollisions among asteroidal and cometary parent bodies 
maintain the observed dust population against losses to 
radiative forces. Since dust production is enhanced by 
gravitational stirring, debris disks systems are natural 
targets for giant planet imaging searches. It has been 
established both theoretically and observationally that 
warm dust is transient. Its presence thus serves as a 
marker for a young or dynamically active planetesimal 
belt, and clearly indicates that the host star possesses 
some kind of planetary system. Like self-luminous exo¬ 
planets, the brightness of a debris disk decays with time. 
To first order, the magnitude of the warm excess serves 
as a chronometer: 22/24 /rm flux densities larger than 
1.2 X the stellar photosphere are almost always found in 
sources younger than 1 Gyr, and more typically, in stars 
with ages less than a few hundred Myr. This theoretical 
expectation has been borne out through b oth observa¬ 
tional work and computational modelling. [Rieke et al.j 
(2005) have shown that warm excess, as observed at 24 
/rm, IS a steeply declining function of stellar age in a sam¬ 
ple of 266 A stars. Models of the dynamical evolution of 
planetesimal swarm s confirm this behaviour for A stars 
([Wyatt et al.||2007D, punctuated by transient spikes of 
higher dust content following major planetesimal colli- 
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TABLE 1 

Properties of HR 3549 


Properties 

value 

Identifiers 

HR 3549, HD 76346, HIP 43620 

Coord, (hms, dms) 

08 53 03.77832 -56 38 58.1462 

Galactic coord, (deg) 

274.35 -07.66 

Spectral type 

AOV 

Distance (pc) 

92.5 ± 2.5 

V mag 

6.01 

L mag 

6.04 ± 0.05 

^la (mas/yr) 

-22.81 ± 0.26 

* cos (5 (mas/yr) 

36.54 ± 0.26 

RV (km/s) 

23.90 ±2.2 

[3.6] - [22]^im 

0.56 ±0.06 

Teff (K) 

10207 ± 347 

log(ff) (cgs) 

4.20 ± 0.14 

nsinz (km/s) 

212 llGlebocki & Gnacinski|2005|| 

Median Bayesian age (Myr) 

230 -- 

68% age Cl (Myr) 

120-360 

95% age Cl (Myr) 

10-390 

2D interpolated age (Myr) 

200"^^^ 

Median Bayesian mass (Mq) 

2.32 

68% mass Cl (Mq) 

2.2-2.5 

95% mass Cl (Mq) 

2.1-2.6 


sions. 

Here we report the detection of a bound low-mass com¬ 
panion to dusty host star H R 3549. This discovery, rem¬ 


iniscent of K Andromedae (jCarson et al. 2013 
et al.|2013| ) 


Hinkley 

adds to the collection of extreme mass ratios 


sub-stellar companions filling the brown dwarf desert. 

2. FUNDAMENTAL PARAMETERS OF HR 3549 

HR 3549 (HD 76346, HIP 43620) is a main sequence 
AO star of visual magnit ude 6 (Table B. T he parallax 
measured by Hipparcos (van Leeuwen P007| is 10.82 ± 
0.27 mas, corresponding to a distance of ~ 92.5 pc. The 
WISE satellite measured W1-W4 = 0.56 ± 0.06 mag of 
excess at 22 /rm, and no excess at 12 ^m (Cutri et al. 


(20131, data from AllWISE cryogenic sky survey), h'or an 
Ao photosphere, this excess corresponds to a minimum 
Ldust/Lstar of 10“^, similar to planet-bearing HR 8799. 
This value is for the case of a blackbody that peaks at 22 
/rm. The presence of any emission extended over a range 
of temperatures would add to this value. 

We used the Bayesian Analysis for Nearby Young As- 
sociatioNs H (BANYAN H) online tool to determine 
the membership p robability of HR 3549 to nearby youn g 
kinematic groups (|Malo et al.||2013| |Gagne et al.||2014D . 
This tool is based on a comparison of Galactic position 
(XYZ) and space velocity (UVW) to well-defined mov¬ 
ing groups closer than 100 pc and younger than 200 Myr. 
No clear association can be found with any known nearby 
young moving group, so HR 3549 is likely a field AO star. 

Our direct interpolation in log (Teff)—log jq) space rela¬ 
tive t o PARSEC vl.l evolutionary models (Bressan et al. 
20121 yields 200 Myr, with a 40-220 Myr 68% confidence 
interval, determined from a distribution of 50,000 inter¬ 
polated ages reached through Monte Carlo propagation 
of the associated errors in log (Tpff) — log (q) . However, 
as pointed out in David & Hillenbrand (2015), direct in¬ 
terpolation does not account for the nonlinear mapping 
of time onto the H-R diagram, nor the non-uniform dis¬ 
tribution of stellar masses observed in the galaxy, and 


can lead to biases. 

3. OBSERVATIONS AND DATA REDUCTION 

We observed HR 3549 with NAOS-CONICA at the 
Very Large Telescope as part of program 090.C-0486(A), 
“L-band adaptive optics imaging of exoplanets around a 
sample of dusty A stars recently discovered by WISE” 
(PI: Mawet) on 2013-01-13. We used the Lp-band fil¬ 
ter. L-band (centered around 3.8 /rm) is a compelling 
and competitive filter for ground-based planet surveys. 
This wavelength range offers significant advantages com¬ 
pared to shorter wavelengths: (i) The L-band contrast 
of planetary-mass companions with respect to their host 
stars is more favorable than in the H and K bands 


(Baraffe et al. 

2003 

) so that lower-mass objects can be 

probed; and (ii 

) longer wavelengths lead to better image 


quality and a more stable point-spread function (PSF), 
with Strehl ratios well above 70%. These advantages 
outweigh the increased sky background in the thermal 
infrared and the loss in resolution. Finally, we note that 
background star contamination probability rates will be 
near-zero at small separations in L-band, so minimal 
follow-up time will be needed to confirm candidates. 

The data were acquired as a sequence of 10 exposures 
on a 5-point dithering pattern with offsets of about 6 
arcseconds each. Each exposure was the average of 100 
frames with 0.25 s integration time, making for a total 
open shutter time of 250 s. We used the pupil track¬ 
ing mode where the instrument co-rotates with the tele¬ 
scope pupil to fix diffraction and speckles to the detec¬ 
tor reference frame, allowing the sky to counter-rotate 


wiLii ine oaranacLic aiiKitJ, eueeiiveiv e 
differential imaging (ADI Marois et al. 

naun 

2006 

ng angular 

). ADI re- 

quires sufficient sky rotation to avoid se 

It-subtraction of 


the companion signal at small angles, often leading to 
long sequences. Our strategy to overcome this limitation 
of ADI has been to limit the sequence duration and rely 
on our uniform target sample to build a library of refer¬ 
ence PSFs, and perform reference star differential imag¬ 
ing (RDI). Indeed, if the sample is uniform in brightness, 
covers a reasonable range of observing conditions and is 
spread uniformly across the sky, correlated speckle pat¬ 
terns can be retrieved in the library of PSFs. 

We reduced the data by subtracting a background 
made out of median-combined dithered frames, dividing 
by a flat field and interpolating for bad pixels and other 
cosmetics. The reduced images were then processed by 
two independent speckle cali bration packages base d on 
principal component analysis (Soummer et al. 2012), and 
the library of reference PSF for the speckle calibration. 

4. DISCOVERY OF A CANDIDATE LOW-MASS 
COMPANION TO HR 3549 

The 2013-01-13 data set showed a point source at a 
separation of ~ 0'.'9, and position angle (PA) of ~ 157°. 
The detection is unambimous with a signal-to-noise ratio 
(SNR) of ~ 6 (see Fig. |H left) The TRILEGAL star- 
count model (Girardi et m. estimates a probability 
of 7 X 10“^ that it is an unrelated background object. 
While this number must be interpreted with caution, it 
is significant enough to warrant follow-up observations. 

We acquired the second epoch data set as part of pro¬ 
gram 094.C-0406(A), “L-band adaptive optics imaging of 
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Fig. 1.— Left: discovery image of HR 3549 B, taken in January 2013, and confirmation image in January 2015. Right: Common proper 
motion analysis of HR3549 B. The black dot shows the relative position of the candidate companion in 2013. The empty circle shows the 
relative position of the point source in 2015 if it were a background object, accounting for parallactic motion (plain curve). The gray dot 
is the measured position of the bound companion in 2015. 


exoplanets around a sample of dusty stars recently dis¬ 
covered by WISE. Part III: candidate follow-ups” (PI: 
Mawet), on 2015-01-14, following a similar strategy as 
for the discovery epoch. The candidate was detected 
at roughly the same location though with a lower SNR 
(~ 3, see Fig. right). The lower SNR appears to be 
mostly due to the characteristics of the Aladdin2 de¬ 
tector installed in CONICA in January 2015 to replace 
the faulty Aladdin3 detector, and the additional thermal 
background induced by the different shielding character¬ 
istics of this detector configuration. 


5. ASTROMETRY AND COMMON PROPER 
MOTION ANALYSIS 

We performed astrometric and photometric measure¬ 
ment of both epochs using forward modelling of the off- 
axis companion point spread function with a Monte- 
Carlo Markov chain (MCMC) sampler in the a, 5, flux 
space. The MCMC sampler model also takes into 
account variable error bars, meaning that systematic 
under-reporting of errors should not affect the final re¬ 
sult. 

The final error bars on the astrometry were then con¬ 
servatively set to the quadratic sum of the MCMC- 
derived error bars, which can still be biased by underly¬ 
ing speckle noise, and the empirically derived infl uence of 
SNR on a, S, flux. Indeed, astrometric precision (Guyon 
et al. 12012 1 is proportional to FWHM/(2 x SNR) in the 
speckle-noise dominated regime, where FWHM is the full 
width at half maximum of the resolution element A/D, 
with observing wavelength A and telescope diameter D. 

As far as the astrometric calibration of CONICA is 
concerned, we used the 27.1 ± 0.04 mas/pixe l plate scale 
and — 0.45 ± 0.09 degree true North offset of |Absil et al. 
(2013) for the 2014 epoch, and the 26.99±0.02 mas/pixel 
plate scale and -1-0.31 ± 0.02 degree true North offset for 
the 2015 epoch, measured after the NACO recommission¬ 
ing with the Aladdin2 chip on UTl (personal communi¬ 
cation from NACO instrument scientist Julien Girard, 
European Southern Observatory). These additional sys¬ 


tematic offsets and corresponding error bars were folded 
into both epoch astrometric positions and associated er¬ 
ror bars (Table ^ to perform the common proper mo¬ 
tion analysis (CPM) shown in Fig. H (right). Following 
the GPM analysis, we determined fnat the probability 
of the discovered candidate to be a background object is 
2 X 10"®. 


6. AGE OF HR 3549 

As mentioned earlier, direct interpolation in an H- 
R diagram or color-magnitude di agram can lead to bi¬ 
ases i n derived ages. Following [David fc Hillenbrand| 
(|2015|) , we thus used a Bayesian analysis of the star’s lo¬ 
cation in log Teff — log g space relative to solar meta llicity 


{Z = 0.015) PA RSEG vl.l evolutionary models (Bres- 
san et al. 2012), yielding a 68% age confidence interval 


ot 120-360 Myr, which is consistent with the average age 
of AO field dwarfs. We performed likelihood calculations 
on a 1000x1000 grid from 1 Myr-10 Gyr in ag e, and 1- 


10 M pt in mass. We used the Salpeter IMF (Salpeter 


1955) as the prior on stellar mass and a uniform prior m 
linear age (i.e. constant star formation rate). The atmo¬ 
spheric parameters are determined from combination of 
uvby0 photome try and ATLAS9 models (jGastelli & Ku- 
2004 2006), and the details of both the atmospheric 


rucz 


characterization and age de termination are presented in 
David & Hillenbrand (2015), hereafter DH15. 

Intermediate-mass stars are rapid rotators, and the ef¬ 
fect of this rotation is to make the star appear cooler, 
more luminous, and hence older than a non-rotating star 
of the same mass. Additionally, rapid rotators spend 
a longer time on the main sequence than slow rotators. 
Recently, DH15 and Brandt & Huang (2015), hereafter 
BH15, have shown that the rotational effects on the in¬ 
ferred ages of intermediate-mass stars are substantial]^ 


® We note these authors account for rotation in different ways 
when determining a ges for such stars. While DH15 apply the cor¬ 
rective formulae of [Figueras Sz Blasi| | |1998| | to the atmospheric 
parameters prior to age determination (assuming vsini ^ Urot), 
BH15 marginalize over projected rotational velocity. 
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T (Myr) 


Fig. 2.— Left: Dete rminations of the star’s location in — logg space. The solid curves are solar-metallicity PARSEC vl.l i sochrones 


[! 


represent 


pectivelyT 


clearly discrepant. 

Right: Posterior probability di stributions in age in bot h the rotation-corrected (red) and uncorrected (black) cases from DH15. The teal 
curve is the age posterior from [Brandt fc Huang] (|2015^. Taken collectively, the age of HR 3549 is in the range of 50-400 Myr. 


Consequently, we explored two possible solutions for 
the age of HR 3549. The first solution is based on 
the rotation-corrected atmospheric parameters derived in 
DH15, Toff = 11176 ±380 K, log^ = 4.33 ±0.14 dex (see 
Figure [^. At this position, direct interpolation yields 
an age of 43^3^ Myr, which is consistent with a previous 
estimate of 55 Myr (Westin 1985P The Bayesian age 
analysis in this case yields a median age of 110 Myr with 
68% and 95% confidence intervals of 10-150, 10-250 Myr, 
respectively. 

The second, more likely, solution is based on the un¬ 
corrected atmospheric parameters from DH15, listed in 
Table [2 This set of atmospheric parameters yields a me¬ 
dian Bayesian age and 68% confidence interval of 230^33° 
Myr, with a corresponding interpolated age of 200^300 
Myr (where the uncertainties are determined from Monte 
Carlo error propagation). We consider this solution more 
likely for the following reasons: (1) the effective tem¬ 
perature is more consistent with moder n spectral type 
scales which suggest Tpff s s 9700 K at AO ( Aclelman|20C)4 


Pecaut & Mamaj ek|2013| ), (2) both Teff and log g in this 


case are consistent with the mean of previous determi¬ 
nations from the literature, (3) this older age is more 
consistent with the average age of AO field dwarfs, (4) as 
noted in § [^ the star is unlikely to be associated with 
known young moving groups, (5) there is some evidence 
presented in DH15 that the procedure used to determine 
the rotation-corrected parameters is over-aggressive, (6) 
there are no other signihcant indicators of youth for 
HR 3549, (7) the older age is c onsistent with other m od- 
ern estimates: 3 00 ± 51 Myr (jZorec fc Royer||2012|) and 
138 ± 98 Myr (Gerbaldi et ^ 1999), both from H-R 


diagram analyses. 

Figure]^ (left) demonstrates the position of HR 3549 in 
log Teff —tog g space, relative to evolutionary models. The 
difficulty of age-dating intermediate-mass stars on and 
near the main sequence is evident from the typically large 


jWestinj ( 1985| discusses age errors in a general sense but does 
not provide intormation for assigning errors to the AO stars in that 
sample. 


uncertainties in surface gravity (or equivalently, luminos¬ 
ity). The effect of including rotation is also illustrated 
by the significantly different atmospheric parameters ob¬ 
tained. 

Figure (right) shows the posterior probability dis¬ 
tribution mnction (PDF) in age for HR 3549, originally 
derived in DH15. Also depicted is a PDF derived from 
the similar Bayesian app roach to isochrone age-dating of 
Brandt & Huang (20151. The BH15 PDF was generated 
assuming a Gaussian prior on [Fe/H], with /r=-0.1 and 
(7=0.2 dex, consistent with the distribution observed for 
intermediate-mass stars in the solar neighborhoodj^ The 
BH15 PDF is broader than that of DH15 due to the fact 
that those authors marginalize over mass, metallicity, ro¬ 
tational velocity, and inclination, three parameters which 
can substantially affect the inferred ages of intermediate- 
mass dwarfs. In contrast, DH15 marginalize over mass 
only, while the typically minor differences in metallicity 
are implicitly accounted for in the atmospheric parame¬ 
ter uncertainties. For comparison, the BH15 PDF yields 
a median age of 220 Myr, with 68% and 95% confidence 
intervals of 10-290, 10-500 Myr, respectively. 

All estimates suggest the age is <500 Myr, and no¬ 
tably the three most recent estimates suggest r > 200 
Myr, which is consistent with the average age of AO field 
dwarfs. The range of published ages illustrates the diffi¬ 
culty of age-dating on and near the main sequence, par¬ 
ticularly for intermediate-mass stars for which empiri¬ 
cal age-dating methods are either non-existent or uncal¬ 
ibrated. We adopt the median Bayesian age and 68% 
confidence interval as the fi nal age and uncertainties fo r 


HR 3549 , T « 230_{3q Myr David & Hillenbrand 
However, the literature range can also be used 


(2015). 
‘to infer 

the loosest reasonable constraints on the system age: 50- 
400 Myr. For completeness, we consider this broad range 
of plausible ages when inferring the companion mass. 

7. PHOTOMETRY AND PROPERTIES OF THE 
COMPANION 


® The BH15 PDF was downloaded 
http: //bayesianstellarpcirameters. info/ 


from 
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TABLE 2 

Astrometry and photometry of the low-mass companion to HR 3549 


Data set 

Filter 

SNR 

Aa (») 

A5 (") 

m 


NACO 2013 

Lp 

~ 6.25 

01'333 ± 07009 

-07806 ± 07009 

13.85 ±0.25 

9.03 ±0.26 

NACO 2015 

Lp 

~ 3.1 

01'334 ± 07015 

-07788 ± 07015 

13.63 ±0.5 

8.5 ±0.505 


The forward modeling MCMC approach to measure 
the companion photometry described above yields an ab¬ 
solute L magnitude of 9.03 ± 0.26 for the companion, 
including parallax/distance and L-band magnitude un¬ 
certainties. To infer the compani on properties , we used 
the BT-Settl evolutionary model (Allard 2014) and our 
adopted conservative age range of 50-400 IVlyr. BT-Settl 
covers the range from solar-mass stars to the latest-type 
T and Y dwarfs, and reproduces the formation of clouds 
and in particular their clearing at the L/T transition. 
We derived a mass range of 15 — 90 Mj (Figure [^, and 
effective temperature between 1900K and 2700K, facing 
HR 3549 B in the L-dwarf regime. 



L mag 


Fig. 3.— Mass (Mj) vs Luminosity (L mag) diagram for the 
evolutionary model BT-Settl. HR3549 B’s L-band photometry and 
associated uncertainty is in shaded red. 


8. DISCUSSION 

No far-infrared observations of HR 3549 are listed in 
the data archives for the Spitzer or Herschel space tele¬ 
scopes, so the properties of its debris disk can only be 
estimated from the WISE survey data. The AllWISE 
survey magnitudes are 6.04±0.05, 5.97±0.04, 6.06±0.02, 
and 5.50 ±0.03 - corresponding to 3.6, 4.5, 11.8, and 22.0 
fxm respectively. An excess of 0.56 ±0.06 mag is detected 
at 22 /im, but no excess is evident at shorter wavelengths. 
Using a 3 tr upper limit to the 11.8 fj,m excess (0.05 mag), 
and assuming blackbody emission, the upper limit to the 
disk dust temperature is 168 K. For grain sizes of a few 
microns, this corresponds to a disk inner edge of r > 10 
AU. 

While the WISE data provide no constraint on the disk 
outer radius, HR 3549 B appears to be massive enough 


to gravitationally clear any disk material in its vicinity. 
It is unlikely that any disk material located beyond the 
companion’s >80 AU orbital separation would be warm 
enough (^ 120 K) to emit significantly at 22 ^m. The 
most likely scenario is therefore that HR 3549 B orbits ex¬ 
terior to a warm dust belt that is the source of the WISE 
excess. An exterior cold dust belt beyond HR 3549 B or¬ 
biting radius may also exist but would require ALMA 
observations to detect. 

The HR 3549 syste m shares several features with the k 
Andr omedae system (Carson et al. 2013 Hinkley et al. 
2013). Like k Andromedae, HR 3549 is a late-B/early- 
A type star with a mass ~3 A/©, and a poorly con¬ 
strained age. Further, the companions in both systems 
are “brown dwarf desert” objects with masses in the 
range of ~15-80 Mj„p, corresponding to mass ratios of 
^1%. Such “extreme mass ratio systems” (e.g. Hinkley 
et al.|201^ ), are particularly important for constraining 
the formation mechanisms of stars in this mass range. 


Indeed, several works (e.g . Delgado-Donate et al. 2004 
Kouwenhoven et aL||2007 ) have suggested that the mul¬ 


tiplicity ot intermediate-mass stars may serve as a reser¬ 
voir for the conserved initial angular momentum in the 
protostellar cloud. 

Lastly, we note that no other companion more massive 
or similar to HR 3549 B is detected around HR 3549 from 
the effective inner working angle of ~ 0”.3 (projected 
separation of ~ 30 AU) to the outer edge of the effective 
field of view of ~ 12” (~ 1000 AU). 


9. CONCLUSION 

This paper presented the detection of a substellar com¬ 
panion orbiting at a projected separation of ~ 80 AU 
around HR 3549, a disk-bearing AOV star. The charac¬ 
terization of the companion is made difficult by the un¬ 
certain age determination for field A stars. Spectroscopic 
follow-up using medium resolution slit spectroscopy is 
the next step to fur ther characterize the low-mass ob¬ 
ject around HR3549 ( Hinkley et al.|201^ ). Indeed spec¬ 
tral indices to quantitatively measure the strength of the 
FeH, VO, KI, spectral features have been demonstrated 
to be robust age markers, but require medium resolution 
spectroscopy (i? > 300). 
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